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ABSTRACT 

Using an X-ray stacking procedure, we provide a robust upper limit to the X-ray luminosity per 
object of a set of 54 z w 5.8 galaxy candidates in the Hubble Ultra Deep Field (HUDF), which is within 
the 1 Ms-exposure Chandra Deep Field-South (CDF-S). With an effective total exposure of 44 Ms for 
the stack, the 3 a flux-density limit of 2.1 x 10~ 17 erg cm -2 s _1 (soft-band) gives a 3 c upper-limit 
luminosity of Lx = 8 x 10 42 ergs _1 per object at a rest-frame hard energy range of ~ 3-14 keV at 
z ~ 5.8 for a photon index of T — 2. For an active accreting black hole (or "mini-quasar" ) emitting at 
the Eddington luminosity (??Edd = 1) ; and the Sazonov et al. average-QSO spectral energy distribution, 
we calculate an upper limit on the black hole mass, M^h. < 3 x 10 6 tj^^Mq (3 a). The X-ray limit 
further implies an upper limit on the rate density of UV ionizing photons from accreting black holes 
at that redshift, ni on i Z < 2 x 10 51 s _1 Mpc~ 3 (3 a), which is less than 1/10 of the number needed to 
ionize the universe. Because the constraint is anchored in the rest-frame hard X-ray regime, a steeper 
spectrum for "mini-quasars" would imply relatively fewer ionizing UV photons. Unless there are 
large populations of active black holes around this mass that are unassociated with luminous galaxies, 
mini-quasars do not appear to contribute significantly to the budget of ionizing photons at z « 6. 
Subject headings: galaxies: formation - galaxies: evolution - galaxies: high-redshift - early universe - 
X-rays: galaxies 



1. INTRODUCTION 

The identification of the epoch or epochs of reioniza- 
tion is not a solved problem, though there is tension 
between the high redshift (z re ion ~ 10-17) of reioniza- 
tion implied by th e WMAP results l|Sr)ergel et alJl2003t 
iKoeut et al.ll2003j) and the lower redshift (z re ion ~ 6) 
suggested by the detection of the Gu nn-Peterson trough 
in z » 6 luminous QS Os by SDSS {Becker et all 1200 It 
iDiorgovski et al.ll2001l) . These observations suggest that 
the transition from an opaque to a fully reionized uni- 
verse must have at least ended by z ss 6. One of the chal- 
lenges faced toda y is to identify the sources of the ioniz- 
ing ph otons (e.g. lLoeb fc Barkanall200H iStiavelli et alJ 
l2004al and references therein), which are still unac- 
counted for from observed stellar and accreting black- 
hole ("active") populations. 

The Advanced Camera 
iBem'tez et all 12004) on the 
(HST) has made routine 
didate galaxies at redshi ft 
120031 iBouwens et al.l 120031 IStanwav et 
Dicki nson et alJ 12004 IGiavalisc o et alJ l2004bl) . These 
objects are identified by an extension of th e "Ly- 
man Break" tec hnique {Steidel fc Hamilton! 119921: 
iSteidel et alj [l996), whereby the opacity from neutral 
hydrogen in the intergalactic medium (or the Lyman-a 
forest) suppresses th e flux at waveleng ths shortward 
of the Lyman break {Madau et alj [l996), or shortward 
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for Surveys (ACS; 
Hubble Space Telescope 
the discov ery of can- 
z pa 6 llBunker et al.l 
t~aiTl2003L 12004 



of the Lyman-a emission {Dickinson et alJ 12004) ). 
The HST/ACS 1775 and z 850 filters bracket Lyman-a 
(A res t = 1216 A) around z ta 5.8, so searches for "i- 
dropouts" will isolate candida t e galaxies around tha t 
redshift {Stanwav et all l2003t Dickinson et~aT] 12004) . 
as well as possible "interlopers" in the guise of cool 
stellar dwarfs or faint evolved galaxies at z w 1-2. 
Some of these high-redshift candidate galaxies have 
been confirmed spectroscopically fr om the ground, and 
from the G800L gr ism of ACS {Pirzkal et all l200l 
iMalhotra et~alll2004) . 

Using data fr om the ACS-based Hu bbl e Ultra Deep 
Field (HUDF; iBeckwith et"all 12001 . iBunker et all 
(2004) identified 54 spatially-resolved i-dropout candi- 
dates to zs5o ~ 28.5. It is still an open question whether 
some combination of these galaxies, possibly with the 
contribution of as yet unveiled populations, have enough 
energetic phot ons to ionize (or keep ionized) the univers e 
at this epoch {Bunker et alJl2004 IStiavelli et al.ll2004bjh 

The need to account for the large number 
of required ultraviolet ph otons for ionization 
{Madau. Haardt. fc ReesI Il999fl has given rise to a 
broad range of proposed sources, ranging from known 
candidates (stellar and active sourc es) to exotic ones 
(e.g. decaying heavy neutrinos, lHansen fc Haimail 
2004). The stellar sources can be divided broadly into 
two categories. T he first-generation "metal-free" stars 
{Abel et alJl2002D . which may arise in M dm > 10 5 M Q 
dark matter halos, are expected to have masses of 
M* £ 60-M© 4 . These stars, though rare, may produce 
enough photons to contribute most or all of the UV 

4 Such stars with masses in the range M* ~ 60— 14OM0, and 
M, > 260M (n are predicted to collapse directly to a black hole 
IHeger et al J 12003) . with M, ~ 140-260 M stars annihilating 
themselves through pair-creation processes. Such collapsed stars 
may be the seeds of the "mini-quasar" black holes discussed below. 
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photons needed to ioniz e the universe even at very 
early times (z ~ 20; e.g. iSomerville fc Livid 120031) . If 
such objects are to be found at z ~ 6, it is unlikely 
that significant amounts of X-ray emission would be 
associated with them. Population II stars, which may 
also contribute to the ionizing radiation budget at high 
redshift, could have associated X-ray emission from 
X-ray bina ries that result from different modes of star 
formation ( Grimm et al.M2003j) . Though such relatively 
locally determined relations may not hold over all time, 
they can be used to explore relevant limits to star 
formation rates at high redshift. 

It is plausible that many of the first collapsed ob- 
jects produced black holes of some mass, that could then 
be sources of energy. This possibility is potentially at- 
tractive as a source of energetic photons, particularly if 
there are mechanisms for self-regulation of their growth 
llWvithe fc Loebl20 03). The mass-increase e-folding time 
of black holes emitting at their Eddington luminosities, 
the Salpeter timescale, is Ts alp rj 1 x 10 8 (e /0.30)- 1 yr. 



TABLE 1 

Results from the X-Ray Stacking 



where e is the energy con version efficiency ijElvis et al l 
l2002UGammie et alJ l2004). The T Salp is independent of 
the black hole mass. If the active lifetime of an accret- 
ing black hole at z ~ 6 is several times Ts a i P , so that 
the probability of any one host galaxy being observed 
while "active" is high, it is possible that the majority of 
galaxies hosting black holes could be observed during an 
active stage. 

In this paper we use one plausible set of dropout- 
selected galaxies at z « 5.8, described in §|5| None of 
these sources are individually detected in the deep Chan- 
dra observations in that field. We therefore construct an 
extremely deep stack of the full sample, as set out in 
§13 and compute an upper limit to the average X-ray 
luminosity for each galaxy. This number is cast against 
proposed sources of ionizing photons at that redshift in 
§0] and the possible utility of these limits are discussed 
in §0 We assume the concordance cosmology param- 
eters of a flat universe with = 1 — Q m = 0.7, and 
H = 70 h 70 kms^Mpc- 1 . 

2. THE HUDF AND HIGH- REDSHIFT CANDIDATES 

The HUDF 5 is a 11.3arcmin 2 ACS-based survey of a 
single extrem ely deep field observed w ithin the GOODS- 
South 6 field ijGiavalisco et al.|l20 04a). The observations 
have been made in the same filters as the GOODS 
project, -B435, VI306, 2775 and z^bo, and reach lOer depths 
of 

2850 ~ 28.5 (AB). This is also the 1 Ms Chandra Deep 
Field-South (CDF-S) field, which reaches on-axis flux 
density limits of 5.5 x 10 -17 ergs -1 cm -2 in the soft-band 
(0.5-2 keV) an d 4.5 x 10~ 16 ergs" 1 cm -2 in the hard- 
band (2-8 keV; iGiacconi et aLll2002]l . 

By virtue of the dropout selection with the «775 and 
Z850 bandpasses, the redshift distribution/selection func- 
tion of the high- redshift galaxies is quite narrow, cen- 
tered around z w 5.8. The candidate list in lBunker et al.l 
(2004) is a lower limit on the total number visible in the 
HUDF, because of the very conservative signal-to- noise 
requirements applied in the detection band. At the same 
time, by the strong color criterion, there should be few 
foreground [z w 2) galaxies. Therefore we adopt the 

5 http://www.stsci.edu/hst/udf/ 

6 http:/ /www. stsci.edu/science/goods/ 
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counts; fractional counts are due to fractional detection cells within 
the apertures; (3) Effective exposure time; (4) Exposure-corrected 
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tabula ted list of candidates directly from iBunker et alJ 
( 2004) , several of which have already been spectroscopi- 
cally confirmed to be at the claimed redshift range. 
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Fig. 1. — The distribution of X-ray counts for all candidate 
2 ~ 5.8 galaxies. The vertical dashed lines mark the mean back- 
ground count levels in each band, as determined from Monte Carlo 
simulations with 10,000 trials. 



3. THE X-RAY STACK PROCEDURE AND RESULTS 

Stacking the X-ray images of individually undetected 
objects has offered exceptional insights into the na- 
ture of star formation and activity in otherwise X-ray- 
invisible extremely distant objects (e.g. [Alexander et alJ 
I2002albt iNandra et alJ l2002: Re ddv fc Steide]ll2004D . In 
the present work, we target the positions of all the 
sources described in the previous section. A search radius 
of 1.5 times the error in the X-ray position was applied 
to match each of the optical positions wit h the X-ray 
source catalogs bv lGiacconi et al.l (|2002(1 and lBauer et alJ 
(2004). For each object position, exposure-corrected 
counts were extracted within the 90% encircled-energy 
(EE) radii in the full (0.5-8 keV), soft (0.5-2 keV), and 
hard (2-8 keV) energy bands. Also at each position, in 
order to quantify the significance (or lack thereof) of the 
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counts measurement, we perform Monte Carlo (MC) sim- 
ulations with 10,000 trials using source-free circular re- 
gions of 20arcsec radius around each galaxy. The net 
aperture counts, then, are obtained by subtracting the 
mean value of the (near-Gaussian) MC distribution from 
the number of average aperture c ounts. Details of th e 
stacking procedure are discussed in llmmler et all $2004). 

Using the above apertures to count the number of X- 
ray photons at each galaxy position, we detect 319 counts 
in total, distributed as shown in Fig. ^ an d tabulated in 
Table ^ The exposure-corrected count rate in the full 
band is ~ 2.6 x 10 _6 s _1 in the soft band, with a total 
effective exposure time of 44.1 Ms. There is no formal 
X-ray detection in the stacked data. Assuming z = 5.8, 
a power law spectral energy distribution (SED) with a 
photon index of T = 2, and a foregro und absorbing co l- 
umn density of N m = 8.0 x 10 19 cm" 2 (|Stark et alll992|) . 
we find a band-corrected 3 a limit on the luminosity per 
galaxy of Lx < 8.3 x 10 42 ergs _1 (3<r) in the rest-frame 
energy band (l + z) x (0.5-2 keV) w 3-14 keV. The con- 
version between counts and flux is not very sensitive to 
the assumed spectral properties — a photon index of 
T = 1.7 instead of T — 2 gives a flux which is ^ 5% 
lower. As the fc-correction for a power law with pho- 
ton index T is k(z) = (1 + z) r ~ 2 , the X-ray luminosity 
constraints can correspondingly be lower by ~ 40% for 
r = 1.7. 




Fig. 2. — The cumulative surface density of X-ray emitting active 
objects versus flux limit. The model lines are drawn for quasar 
number co unts at redshifts of z = 5 (dashed line) and 2 = 7 (solid 
line) from Wvithc & Loeb (2003). The limits determined from 
this study for z 5.8 are shown (black box). The counts fr om 
IBarger et all 12003) at z > 5 (diamond) . ICristiani et alj 120041) at 
3.5 ^ z 5.2 (triangle), and lFan et alj |200^) at z }t 5.7 (star) are 
shown for comparison. 

4. HIGH-REDSHIFT SOURCES OF IONIZING PHOTONS 

There is still no consensus on the source or sources 
of ionizing photons at redshifts z ~ 6 and beyond (see 
lLoeb fc Barkanal 120011 for a review of candidate pop- 
ulations), thoug h it is po s sible that a complex reion- 
ization history ijCenl 120031 IHui fc Haiirianl I2003j) could 
have been driven by different combinations of ionizing 
sources at different t i mes (e . g. ISo mervill e fc Livioll2003t 
ISokasian et"aD l200l IMadau et alJ I2004|) . The 

highest redshift where there are presently extensive data 
is z w 6. The success of identifying how the universe 
stays ionized at that redshift is mixed ijYan fc Windhorst! 



120041 IBunker et al.ll200l IStiavelli et al.M2004bt) . In this 
section, we consider the implications of our measured 
limits. 

4.1. Star formation 

Without makin g any corrections for obscuration, the 
iKennicuttJ l)1998j) relation SFR UV (M yr" 1 ) = 1.4 x 
10~ 28 i,y(ergs~ 1 Hz -1 ) (for a Salpeter IMF) gives a typ- 
ical SFRuv ~ 10A^©yr _1 for each galaxy, where we 
have used the z 850 photometry (A rcst ~ 1400 A). At 
low redshift, there is a well-explored association between 
star formation rates and X-ray luminosity, driven by 
the aggregate effect of low- and high-mass X-ray bi- 
naries (LMXBs and HMXBs, respectively), which are 
formed at different stages of stellar evolution. At 
star formation rates above se veral Mfp yr" 1 , promptly- 
formed HMXBs dominate, and iGrimm etaLl 1120031) find 
SFR X (M yr- 1 ) = L 2 ^ 8keV / (6.7 x 10 39 ergs" 1 ). Using 
the rest-frame hard X-ray limit, we calculate SFRx < 
llOOMQyr -1 (3cr). This limit is therefore not use- 
ful compared to the measurements fro m the rest-UV 
(Bun ker et al .120041 IStiavelli et all200 4b). and we do not 
discuss it further. 

4.2. Accreting black holes 

There is theoretical motivation for black hole 
accretion-powered "mini-quasars" at high redshift (e.g. 
IMadau et all 12004). These could begin for ming at very 
early times (z ^ 20) from the collapse of first-generation 
massive stars, and grow substantially in mass by z sa 
6. Whereas the comoving space density of supermas- 
sive (Mbh ^ lO 9 ^^) black holes at this redshift is ex- 
ceedingly small ijFan et al.ll2003j) . there are only infer- 
ences on the numbe r s of lower-mass bla ck holes (e.g. 
IWvithe fc Lo^l27)0!lBromlev et alJl200l . 

With this noted, however, we are encouraged to com- 
pute what black hole masses the X-ray limits in Tabled 
imply, assuming that these putative black holes are ac- 
tually active. For the purposes of this calculation, we 
scale the bolomctric luminosity to the Eddington lumi- 
nosity, as Lboi = ?7Edd^Edd- Since our observations are 
in a r estricted energy band , we use a fiducial SED tem- 
plate (|Sazonov et alJ 12004]) to calculate the fraction of 
the total energy in the relevant band, ft, (the inverse of 
which is the usual bolometric correction factor). Then 
the luminosity in band b is Lx,b = /b VEdd ^Edd, or 
L x , b = 2.52 x 10 42 (/b/0.02) (^Edd/T-O) x 

(M bh /10 6 A/ Q )ergs- 1 , (1) 
where /b ~ 0.02 is calculated for the observed soft-band 
at z ~ 5.8 directly from t he template . Assuming the best 
case of ??Edd ~ 1.0 (e.g. iFlovdl 12001). we find an upper 
limit on the average mass of active black holes in each 
galaxy of 

M bh < 3.3 x 10 6 (/b/0.02)- 1 (Tjte dd/1.0)- 1 M ® ( 3ff )- ( 2 ) 
Integrating the lSazonov et all l)2004|) template between 
1-4 ryd (between the Hi and Hen edges), we find an 
energy fraction of /uv ~ 0.056. The production rate 
limit of H i-ionizing photons per galaxy is then -/Vj on i z < 
ix (/uv//b)/(M ~ 4.0 x 10 53 s -1 . Considering the en- 
semble of galaxies and the approximate relevant volume, 
we estimate the ionizing photon rate density 

n ioniz < 2.0 x 10 51 s" 1 Mpc~ 3 (3 a). (3) 
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5. DISCUSSION 

If remnants from z ss 20 M s t ~ 100 M© population III 
stars form the seeds of early black holes, which then 
grow through radiative accretion until z « 6, they would 
have masses Mbh ~ 10 5_6 M©, even for a large e w 0.3 
and before considering merging. However, at the limit 
of Eq. |21 the density of (active) black hole mass in- 
ferred is pbh.6 < 1-5 x 10 4 M© Mpc~ 3 , which is much 
smaller than the present-day density of all black holes, 
p hh = 2.9 x 10 5 M©Mpc" 3 (|Yu fc Tremaind ISoffl . 
Therefore, the present limit is plausible, while allowing 
great latitude for black holes to continue to grow (and 
shine) . 

For full (hydrogen) ionization, the ionization rate 
(Eq. 13 must exceed the recombination rate, n reC omb = 
C ae nn(z) 2 7 where C = (rig}/(n e ) 2 is the dumpiness of 
the ionized gas, qb ~ 2.6 x 10 _13 cm 3 s _1 is the (case B) 
recombination coefficient for a temperature T rj 10 4 K, 
and nn(z) is the particle density of hydrogen. Then, 

Urecorob — C Q^B [(1 - Ihc) Tlb.O (1 + ^) 3 ] 2 

= 3.2 x 10 52 (C/1.0)s- 1 Mpc" 3 , (4) 

where n b ,o = 2.7 x lO^cm -3 ijSpergel et all 12003ft is 
the present-day density of baryons, and Yh c ~ 0.24 is 
the primordial He fraction. Estimates an d simulations 
suggest that the dumpiness is C « 1-30 ijMadau et al.l 
1999; iCenl I2003D . Therefore, even at the upper limit 
of Eq. |3 active black holes associated with luminous 
galaxies cannot be responsible for but a small fraction 
of the total ionizing flux at z » 6. Most changes to 
the calculations would be in the sense of decreasing the 
/ti nt1 y/n rprntT ih ra ti°- Cosmi c variance ijBarkana fc Loebl 
120031 [Somervil le et alJ l2004) could affect these results by 
a factor of two. The surface density of the z « 6 galax- 
ies used here is, h owever, consistent w ith results from 
independent fields ijBouwens et alJl20 03) . 

It is worth noting that if the SED of an ac- 
creting Mhh ,$ 10 M w black hole is harder than the 
iSazonov et al.l fl2004T) temp late, as has been argued 
e.g. bv IMadaJTeHn $M)M, the X -ray limit derived 
above would imply fewer ionizing UV photons. Such 



sources may still be playing an important role with par- 
tial ionization (especially if they are not explicitly asso- 
ciated with luminous galaxies), or may be instrumental 
in a chieving the early reionization suggested b y WMAP 
(e.g. lSomerville et alJ feOOS: Ma dau et al]l2004D 

Our results are consistent with iDiikstra et al.l {2004), 
who use the observed soft X-ray background to infer 
the maximum contribution to the rest-frame hard X-ray 
background at z > 6, and based on strong observational 
limits, rule out a significant contribution to the ionizing 
flux by active accreting black holes. 

Based on a self-regulated prescriptio n for the growth 
of black holes in a ACD M framework. IWvithe fc Loebl 
( 2003) motivate the local iMagorrian et alj l)1998f) rela- 
tion, and make predictions for massive black holes as 
a function of redshift. Their predicted number counts 
of active black holes at high redshift are shown in Fig- 
ure |2 with the data from this work. This shows that 
useful constraints on such ACDM models will be pos- 
sible with only a few times larger survey areas, even 
at sensitivities compar able to the present CDF-S. The 
IWvithe fc Loebl 1I200I black hole mass function (see 
their fig. 5) predicts a cumulative black hole number den- 
sity of n bh ~ 9 x 10~ 4 Mpc -3 down to M hh « 3 x 10 6 M©, 
which is about five times smaller than the upper limit 
from our sample. This suggests that, distinct from the 
census of ionizing photons from such sources, larger deep 
Chandra surveys in fields with large high-redshift galaxy 
samples may achieve stronger constraints on the mod- 
els. Future X-ray missions, such as Constellation- X and 
XEUS will easily reach the relevant sensitivities - if any 
black holes driving "mini-quasars" at z « 6 are there to 
be found at all. 
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